Density dependent electrical conductivity in suspended graphene: Approaching the 

Dirac point in transport 
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We theoretically consider, comparing with the existing experimental literature, the electrical 
conductivity of gated monolayer graphene as a function of carrier density, temperature, and disorder 
in order to assess the prospects of accessing the Dirac point using transport studies in high-quality 
suspended graphene. We show that the temperature dependence of graphene conductivity around 
the charge neutrality point provides information about how close the system can approach the Dirac 
point although competition between long-range and short-range disorder as well as between diffusive 
and ballistic transport may considerably complicate the picture. We also find that acoustic phonon 
scattering contribution to the graphene resistivity is always relevant at the Dirac point in contrast 
to higher density situations where the acoustic phonon contribution to the resistivity is strongly 
suppressed at the low temperature Bloch-Griineisen regime. We provide detailed numerical results 
for temperature and density dependent conductivity for suspended graphene. 



I. INTRODUCTION 

It has been shownir— that suspended graphene can 
achieve very high mobilities since various annealing tech- 
niques can remove much of the extrinsic impurities un- 
avoidably present in graphene on substrates^—. Such 
ultrapure graphene (in this article, 'graphene' would 
mostly imply 'suspended graphene' without any sub- 
strates) with ultrahigh mobility is of considerable im- 
portance for a number of reasons. First, a careful 
comparison between graphene experimental data with 
and without substrates could inform the community 
a great deal about the type of disorder operational 
in graphene on various substrates and the associated 
scattering mechanisms limiting graphene mobility on 
substrates 1 ^, thus helping the eventual technological ap- 
plication of graphene based devices. Second, ultrapure 
suspended graphene enables the study of interaction 
effects&ii~— without the considerable complications aris- 
ing from disorder, and additionally, the absence of dielec- 
tric screening by the substrate enhances the Coulomb 
interaction. Third, as a direct consequence of the high 
sample purity, suspended graphene is a convenient sys- 
tem for the experimental study of the fractional quantum 
Hall effect^. Fourth, the (relative) absence of disor- 
der in suspended graphene (SG) suppresses the electron- 
hole puddle formation around the charge-neutrality point 
(CNP), i.e., the Dirac point, making it relatively easier 
to access the intrinsic Dirac point physics by lowering the 
carrier density^. 

The last item ( "accessing the Dirac point" ) along with 
the first item ("understanding transport in SG") provide 
the motivation of our theoretical research presented in 
the current paper, where we carry out a detailed quan- 
titative study of SG carrier transport as a systematic 
function of temperature and carrier density neglecting all 
complications arising from the inhomogencous electron- 



hole puddles^— in the system (which typically become 
operational below a typical carrier density n c ~ 10 12 
cm~ 2 for graphene on Si02 substrate a 14 ' 15 ). Our the- 
oretical results provide a direct estimate of the disorder- 
limited SG transport properties down to low carrier den- 
sities which experiments should be able to access in clean 
SG samples where inhomogeneous puddle formation is 
pushed down to very low carrier densities. Our work, 
therefore, should provide a benchmark for understanding 
SG transport data as well as for figuring out how close to 
the Dirac point specific SG experimental samples manage 
to approach. 

A characteristic and universal feature of graphene 
transport is the minimum conductivity phenomenon 
where at some disorder-dependent low carrier density 
(n c ) the conductivity shows an approximate saturation 
as a function of carrier density, forming a rough min- 
imum conductivity plateau around the Dirac pointer— 
with a characteristic electron- hole density width of ±n c . 
The characteristic density cut-off n c defining this mini- 
mum conductivity plateau roughly defines how close in 
density [or in energy e c w ep(ji c ) = hvpyjirric, where 
e c is the graphene Fermi energy for carrier density n c , 
and vf is the graphene Fermi velocity] the particular 
graphene sample approaches the Dirac point. Larger the 
n c (or Ec), further the system is from the Dirac point 
no matter how low one tunes the gate voltage since the 
Dirac point is defined only to the uncertainty of n c . It 
is now reasonably well-establishcd ^ 10 ' 17 i 18 that the min- 
imum conductivity plateau and the characteristic den- 
sity cut-off arises from disorder-induced density inhomo- 
geneity (or equivalcntly electron-hole puddle formation) 
in the system which makes it impossible to access the 
Dirac point nominally existing precisely at zero carrier 
density. Instead the disorder-induced density fluctua- 
tions characterized by n c make the zero-density (and as 
such, measure-zero) Dirac point ill-defined over a scale 
of n c . Smaller n c is, closer one can approach the Dirac 
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point by tuning the gate voltage induced carrier density. 
Thus, transport measurements, which probe n c directly 
by definition (since for |n| < n c the conductivity approx- 
imately saturates) provide a clear signature for how close 
to the Dirac point one is able to approach in a particular 
graphene sample. In high-quality SG, n c ~ 10 8 cm~ 2 
can be achieved^, indicating that experiments can assess 
the Dirac point within e c ~ OAmeV ~ 5K. With further 
improvement in SG sample quality, it is conceivable that 
the SG Dirac point could be accessed within 0.5K leading 
to the possibility of studying intrinsic interaction phe- 
nomenon associated with the non-Fermi liquid aspects of 
the Dirac pointier— . This Dirac point accessibility is the 
primary motivation for our detailed current study of SG 
transport properties as functions of carrier density and 
temperature. In this work we assume n c = 0, and our 
results therefore only apply to ultrapure SG samples at 
doping densities above the conductivity minima and/or 
at ksT > e c . 

In addition to the Dirac point accessibility issue dis- 
cussed above, a secondary motivation of our work is 
a qualitative theoretical understanding of realistic SG 
transport in order to assess whether the current exper- 
imental SG samples are in the ballistic or the diffusive 
regime. Several recent SG experimental investigations^— 
conclude that their studied samples are in the ballistic 
regime based on the estimated transport mean free path 
being longer than (or comparable to) the linear sample 
size. Such very long mean free paths imply essentially no 
carrier scattering within the sample (and consequently, 
almost no disorder), and thus the issue of diffusive versus 
ballistic transport in SG samples is an important topic 
of considerable interest to the community. We find that 
this is also a very subtle topic since the extraction of the 
mean free path from the measured conductivity is quite 
nontrivial at low gate voltage (i.e., near the Dirac point) 
where intrinsic thermal carrier occupancy (because of the 
zero band gap nature of graphene) effects become cru- 
cial, and a naive estimate of the mean free path using 
simply the gate-induced carrier density would seriously 
over-estimate the mean free path. In fact, we believe that 
at low carrier densities it is much more sensible to discuss 
the physics simply in terms of the dimensionlcss 2D con- 
ductivity (in units of e 2 /h) rather than in terms of the 
mean free path and/or mobility which are both derived 
by dividing the measured conductivity by a putative car- 
rier density subject to large errors near the CNP. We 
find that a theoretical description based on purely diffu- 
sive transport using the semi-classical Drude-Boltzmann 
theory gives a reasonable description for the experimen- 
tally observed SG transport properties. We believe that 
the only way to definitively establish ballistic SG trans- 
port is to experimentally observe the explicit sample size 
dependent conductivity characterizing ballistic transport 
where conductance and not conductivity is the meaning- 
ful physical quantity, which, to the best of our knowledge, 
has not yet been seen in any SG samples by any exper- 
imental group. We therefore contend, based on our the- 



oretical results, that the currently existing SG samples 
are all high-mobility diffusive samples. 

We consider primarily disorder-induced resistive scat- 
tering in our thcor y 10 i 17 i 18 since our interest is mainly 
the issue of approaching the Dirac point in high-quality 
SG. The phonon effects have been considered elsewhere in 
details^—, and it is straightforward to include phonons 
in the theory, and we do provide some results including 
phonon scattering in the theory since their effect could be 
important at higher (lower) temperatures (carrier densi- 
ties). 

The rest of this article is organized as follows. In sec- 
tion II we describe our basic transport model and pro- 
vide the expected theoretical results for finite tempera- 
ture Drudc transport of intrinsic (i.e. undoped) graphene 
precisely at the Dirac point, which serves as the starting 
point for later discussions. In section III we provide our 
full theory, and then in section IV we provide our nu- 
merical results, concluding in section V with a discussion 
and a summary. 



II. INTRINSIC TRANSPORT AT THE DIRAC 
POINT 

Precisely at the Dirac point (n = 0), assuming no 
disorder-induced electron-hole puddles and T = 0, it 
is easy to see that the semiclassical Drude-Boltzmann 
conductivity ctd at the Dirac point (or equivalently, the 
charge neutrality point) is precisely zero (i.e. infinite 
resistivity) because of the trivial reason that there are 
no carriers to carry any current. We note that in our 
zcroth-order Drude-Boltzmann transport theory the ma- 
trix element of the off-diagonal terms vanishes due to the 
conservation of energy, which gives rise to the zero con- 
ductivity at T=0 and n=0. But more rigorous transport 
theories (such as Kubo formula or self-consistent Boltz- 
mann transport theory which are beyond the Boltzmann 
theory) produce the non-vanishing matrix elements be- 
tween off-diagonal terms even at n=0. Thus the well 
known minimum conductivity appears in these theories. 
Since our analysis is totally based on the zeroth order 
Boltzmann theory the conductivity vanishes for n=0 even 
for chiral graphene. This trivial result is unstable because 
there will be a finite conductivity the moment the car- 
rier density deviates from the precise measure-zero n = 
constraint which is bound to happen at T ^ even at the 
Dirac point by virtue of the low-energy thermal electron- 
hole excitations capable of carrying the current. Unlike 
ordinary band insulators with finite band gaps 2 ^, there 
is no exponential suppression of finite-temperature band 
conductivity in graphene because of its gaplessness. In- 
stead, as is well-known and discussed in some details be- 
low, (Jd(T) at the Dirac point of graphene manifests a 
power law "insulating" temperature dependence, which 
should distinguish the Dirac point behavior from the sat- 
urated conductivity behavior in the presence of electron- 
hole puddle o 18 i 24 . This power-law "insulating" behavior 
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associated with the Dirac point has nothing to do with 
Anderson localization physics and arises entirely within 
the metallic Drude-Boltzmann diffusive transport theory 
because graphene is a gapless semiconductor. 

Consider undoped graphene in the absence of disorder 
(or electron hole puddles), i.e., the chemical potential at 
T = lies at the Dirac point. Then the thermally excited 
number of electrons (and holes) at finite temperatures 
can be calculated 



n = / D(e)np(e)de, 



(1) 



where np is the Fermi distribution function and D(e) = 
<7£/27T7 2 is the density of states of graphene with the total 
degeneracy g = 4 arising from spin (2) and valley (2)and 
7 = Hvf- The induced carrier density at T becomes 



9 ^ (k B Tf 



2tt 12 



r 



T 2 x 0.89 x 10 6 cm" 2 . (2) 



At T — 3QQK we have n = 8 x 10 10 cm" 2 . Thus if the 
conductivity is simply proportional to the carrier density, 
it increases quadratically with temperature. We show be- 
low that the actual temperature dependence of conduc- 
tivity depends on the scattering mechanism. 

In the presence of disorder induced momentum scatter- 
ing the conductivity can be calculated within Boltzmann 
transport theory. In this theory the puddle effect is not 
considered i.e., the theory is valid only for n > n c . The 
conductivity is given 



deD{e)r{e) 



de 



(3) 



where r is the disorder induced transport scattering time. 
Note that in this equation the conductivity is not related 
explicitly to the carrier density. If we assume a constant 
scattering time, i.e. no energy and temperature depen- 
dence of the scattering time, then we have, 



e 2 2 In 2 

h 7 



(T v F )(k B T), 



(4) 



where r is the constant scattering time and the mean free 
path is given by I = r a vp. In this case the conductivity 
increases linearly with temperature. 

Now consider a generalized scattering time. Within 
the Fermi golden rule we have 

- = IT n i E l^( fe < fc ')| 2 (1 - cob - e k ,), (5) 
k> 

where rii is the impurity concentration and Vi is the 
carrier-impurity scattering potential. 

For a short range potential (i.e. 5-range potential) with 
the strength Vi = V$ and the impurity density jij = ns 
we have 



1 _ n s Vg v F £ 
r(e) 4 7 3 



(6) 



With Eq. (J5J we have 



4e 2 7 2 
h n s V 2 



(7) 



Thus the conductivity is independent of the temperature 
for <5-correlated zero-range disorder. 

For unscreened long ranged Coulomb potential, 
Vi(q) = 27re 2 /ft<7 where k is the background lattice di- 
electric constant (taken to be unity for SG), we have 
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where r s = e 2 / is the graphene fine structure coupling 
constant. With Eq. Q we have 



e 2 1 1 



■{k B Tf 



h 3rij r 2 7 2 

For screened long range Coulomb potential, i.e., 

vf ^ 27re2 
Vi(q) 



nqe(q) ' 

where the dielectric function, e(q), is given by 

2-7TP 2 

e(q) = l+f^U(q,T), 
nq 



(9) 



(10) 



(11) 



where IL(q) is the polarizability depending on the wave 
vector and temperature 2 ^. Within RPA we have, 



n(?,T ) 



4^ 



7T/37 2 



In 2 



^ 2 v/RV^F 



dy 



(12) 



where ft = 1/fcgT and e q = hvFq- Finally, the conduc- 
tivity can be calculated to be (with some straightforward 
algebra) 



e 2 1 (k B T) 2 Tl 

a D{T) = —— 272- J ( r 0)> 



h 2irn 



where 



1 + 7rr s /2' 



(13) 



(14) 



and I(tq) is a function which is independent of the tem- 
perature and given by 



7M = y o dtt 2 r(t,r o)vTW 



where 



r(t,r ) 



dx 



e (2tx,r Q y 



(15) 



(16) 



and 



e (z,r Q ) = H 



ln2- 



z / 2 VI - (2gg 
1 + &v 



dy 
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(17) 

Thus the Dirac point conductivity increases quadrati- 
cally with temperature for screened Coulomb potential 
disorder similar to the bare Coulomb disorder results in 
Eq. ©. 

Finally, for scattering of the thermally excited carri- 
ers by the deformation potential coupling to the acoustic 
phonons^ we get the following expression in the high- 
temperature nondegenerate equipartition phonon distri- 
bution regime: 



O'ph 



e 2 8p m vl h 7 2 1 
h D 2 knT 1 



(18) 



where D is the deformation potential, p m the graphene 
mass density, and v p h the phonon velocity. Thus, the 
conductivity decreases inverse linearly with increasing 
temperature. For low temperatures, T < Tbg where 
the Tbg is the so-called Bloch-Griineisen temperature, 
phonon scattering is very strongly suppressed^ and is 
not of any interest in the current work. 

We note that, as expected, the above Boltzmann the- 
oretical semiclassical description of the Dirac point con- 
ductivity, which neglects all interactions 2 ^ and interfer- 
ence effects 2 ^— (but includes thermal excitation, screen- 
ing, and scattering effects quantum- mechanically) , gives 
<jd{T = 0) = at the Dirac point, and the finite <Jd{T) 
for T 7^ arises entirely from the finite density of the 
thermal electron- hole excitations [c.f., Eq. (JTJ] in gap- 
less graphene. The temperature dependence of the finite- 
temperature Dirac point conductivity is entirely a power 
law with ctd(T) ~ T a where a = 0, 1, 2 respectively de- 
pending on whether the scattering mechanism is short- 
ranged or energy-independent or long-ranged (including 
screened Coulomb scattering). In addition, a = — 1 for 
phonon scattering as shown in Eq. (|18p , and in the pres- 
ence of all possible scattering mechanisms, the actual 
temperature dependent Dirac point conductivity would 
be nonuniversal and complex, depending on the strength 
of the various scattering processes in the particular sam- 
ple. It is then easy to see that the experimentally mea- 
sured a exponent could be any number between and 
2, depending on the manifestly non- universal strength of 
various scattering mechanisms in the system. 

The important point to note is that the temperature 
dependence is never exponential, a key qualitative fea- 
ture which helps to distinguish the Dirac point ther- 
mally induced conductivity from the Anderson (strong) 
localization (or gap-induced insulating) behavior. This 
power-law temperature-dependence feature remains valid 
even in the presence of phonon scattering which leads 
to a metallic conductivity (i.e., a negative a) with 
a temperature-dependent conductivity with a power 
law between 1 and 5^ depending on whether Bloch- 
Griineisen regime is relevant or not. We emphasize that 
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FIG. 1. (color online) The Dirac point conductivity gd as 
a function of temperature. <7j (as) indicates the conductiv- 
ity due to screened charged Coulomb disorder with an impu- 
rity density n* = 0.3 x 10 10 cm -2 [short-range disorder with 
Vi = l.5(eVA) 2 ]. a p h represents the conductivity limited 
by acoustic phonon scattering. a to t is the total conductivity 
including all scattering mechanisms. The dashed line indi- 
cates the conductivity due to the Coulomb disorder and the 
short-range disorder. 



phonons can only induce metallic behavior (with the con- 
ductivity decreasing with increasing temperature), and 
as such disorder and phonon scattering together may 
produce a complicated non-monotonic temperature de- 
pendence. In Fig. Q] we show the calculated <jd includ- 
ing all three scattering mechanisms (i.e. short-range and 
long-range disorder as well as acoustic phonons). The 
total conductivity (atot) shows the non-monotonic tem- 
perature dependence. As temperature increases, cr tot in- 
creases due to the dominance of Coulomb disorder at 
low temperatures, but after reaching a maximum con- 
ductivity it decreases with increasing temperature due 
to phonon scattering. This crossover temperature scale 
for up (T) depends sensitively on the amount of Coulomb 
disorder in the system, and increases (decreases) with in- 
creasing (decreasing) Coulomb disorder. We note that 
if Coulomb disorder is weak or absent, &d{T) decreases 
monotonically with increasing temperature because of 
phonon scattering. The fact that experimental low- 
density SG transport datai~— show a nonmonotonic tem- 
perature dependence in the low-density SG conductiv- 
ity clearly indicates that Coulomb disorder dominates 
even the currently existing SG samples (and not just the 
graphene on substrate samples). 

A key aspect of the temperature-dependent Dirac point 
electrical conductivity derived above, which, although 
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rather obvious, has not been much discussed in the lit- 
erature, is that the intrinsic Dirac point behavior is re- 
ally a high-temperature phenomenon rather than a low- 
temperature one since one must have n(T) > n c in or- 
der to see the intrinsic behavior (where n c is the charac- 
teristic cut-off density defining electron-hole puddle for- 
mation in the system). Thus, the intrinsic Dirac point 
physics can only be accessed for T ^s> T c ss 10~ 3 */n^ 
with n c (T) measured in units of cm -2 (K), and the in- 
trinsic Dirac point behavior is completely suppressed by 
the extrinsic inhomogeneous carrier density fluctuations 
associated with the electron-hole puddles. For the ex- 
tremely low value of n c ~ 10 8 cm -2 , we get T c = 10K 
whereas for the usual graphene on SiC>2 substrates, where 
n c ps 10 12 cm" 2 , T c ps lOOOif ! Thus, the intrinsic Dirac 
point conductivity (and its strongly insulating temper- 
ature dependence arising from Coulomb disorder) can 
never be observed in most graphene on SiC>2 samples 
studied in most laboratories, and indeed, in spite of clear 
theoretical predictions for the insulating temperature de- 
pendence of the low-density graphene conductivity 2 ^, for 
a long time it was believed that graphene conductivity 
is essentially temperature independent upto room tem- 
perature (since the electron-phonon coupling constant is 
small in graphene, even phonon-induce metallic tempera- 
ture dependence is fairly weak in graphene at high carrier 
density). 

To observe the intrinsic Dirac point physics ctd(T) 
at low temperatures (< 100 mK) so that various pre- 
dicted interaction induced Dirac point reconstruction 
(or instability)^ 3 - can be experimentally observed (since 
higher temperature strongly suppresses interaction ef- 
fects), one would have to produce SG samples of rather 
extraordinary purity with the puddle induced density in- 
homogeneity being less than 10 4 cm -2 . This seems to be 
a rather daunting task, and it is therefore safe to say that 
the T — > intrinsic Dirac point conductivity is unlikely 
to be experimentally explored in the near future, making 
our current work, where wc consider finitc-tcmperaturc 
Dirac point transport neglecting interaction effects, to be 
relevant for all experimental Dirac point transport stud- 
ies in the near future. 

Before concluding this section we emphasize that we 
are only considering T disorder-limited Boltzmann 
conductivity in our theory, neglecting all interaction ef- 
fects, and for T = our Dirac point conductivity is 
trivially zero. A completely different approach is nec- 
essary to discuss intrinsic Dirac point conductivity in 
clean graphene at T = where interaction and quantum 
interference effects would be important. Such a theory 
is beyond the scope of our work and is not of interest 
to us since we know of no experimental relevance of the 
T = Dirac point conductivity. Second, the inclusion of 
phonon effects is extremely important for the Dirac point 
conductivity behavior as a function of temperature since 
phonons lead to the non-monotonic ctd(T) with metallic 
behavior (dap/dT < 0) at higher temperatures replacing 
the insulating behavior (dao/dT > 0) at lower temper- 



atures with a nonunivcrsal disorder-dependent crossover 
behavior. 



III. CONDUCTIVITY OF SUSPENDED 
GRAPHENE 

Our theoretical model assumes the absence of pud- 
dles in the system, and as such, the theory is cut off 
at some samples-dependent characteristic carrier density 
n c below which the inhomogeneous density fluctuations 
around the Dirac point become important, leading to 
an observable minimum conductivity plateau formation. 
The theory is valid only for n > n c or in the non-plateau 
regime by definition. If n c is very small, as has been 
claimed in several recent experimental studies*—, then 
our theory would apply down to very low carrier density 
(as long as interference and interaction effects are negli- 
gible). Puddle effects on graphene transport properties 
have been theoretically studied elsewhere 2 ^, and puddles 
would introduce additional nontrivial temperature and 
density dependence for n < n c , which is not of interest 
to us in the current work where n c is very small by virtue 
of the ultra-clean nature of SG in general. 

In our theory, we consider three distinct scattering 
mechanisms contributing to the SG resistivity. These 
arc charged impurity, short-range disorder, and in- 
plane acoustic phonon scattering processes. There can 
be other types of scattering mechanisms contributing 
to the graphene resistivity such as resonant scatter- 
ing center a 31 ' 32 , ripples 3 - 3 -, and flexural phonons 3 ^—. 
The short-ranged disorders considered in the resonant 
scatterer o 31 ' 32 modify the density of states of graphene 
(i.e., there are a resonant DOS peak at the Dirac point 
due to the disorders). In this case the short-ranged disor- 
der scattering gives rise to the density dependent conduc- 
tivity. In our Boltzmann theory, to keep the theory to be 
consistent for all other disorders we used the bare DOS of 
pure graphene and we have the density independent con- 
ductivity from the short-range disorder. In this paper 
the short-range disorder represents the zero-range disor- 
der, i.e., V(r) — VsS(r). The calculated scattering time 
with the finite width potential (i.e., V(r) = Vg6(r — ro)) 
does not significantly modify that with the zero-range 
potential as long as r < 2a where a is the lattice con- 
stant of graphene. In suspended graphene the flexural 
phonons may dominate the phonon contribution to the 
resistivit y 35 ' 36 . However, in suspended graphene under 
specific tension induced by contacts (this is the case for 
all available suspended graphene samples) the flexural 
phonon contribution to the conductivity is severely sup- 
pressed and as a consequence the in-plane phonon is the 
dominant scattering mechanism^ 6 -. Since we consider sus- 
pended graphene under the tension our calculated results 
is not affected by the flexural phonons. The possibility of 
still other unknown scattering mechanisms (such as scat- 
tering from the hybridization of electron-hole excitations 
and out-of-plane optical phonons 3 -^) contributing to the 
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graphenc resistivity cannot be ruled out either. But we 
neglect these scattering mechanisms because we want to 
keep the number of parameters to a minimum and also 
because the very high SG mobility and quality imply that 
the overall scattering contributions are small. 

The Drudc-Boltzmann conductivity theory for extrin- 
sic graphene in the presence of induced carriers is a 
straightforward generalization of the theory provided in 
Sec. II except the total carrier density now has an exter- 
nally tunable (through the gate voltage) density in addi- 
tion to the thermally excited intrinsic carriers considered 
in Sec. II. This theory has been much discussed in the 
literature, and we provide below the working equations 
for different contributions to the SG resistivity from the 
three different scattering mechanisms considered in our 
work. The theory below is a straightforward generaliza- 
tion of the theory for the finite-temperature Dirac point 
conductivity developed in the last section. 

Within the Boltzmann transport theor y 10 ' 25 , the con- 
ductivity <j(n, T) is given within the relaxation time ap- 
proximation by 



deD(e)v 2 r(e) 



df{e) 
de 



(19) 



where /(e) is the relevant distribution function. The re- 
laxation time r(e) = r(efe) is given after ensemble aver- 
aging over random disorder configuration by 

1 2tt^ f „i f d 2 k' 

Ji 



-W(e fe ) 



(a) 



(2tt) 2 



x (l-cos0 w )<5(£ fe -4): (20) 
where 9kk' is the scattering angle and is the poten- 
tial disorder causing the scattering with n\ a ^ being the 
2D density of the random impurities (or defects) produc- 
ing the disorder and '( a )' being a label indicating the 
kind of scatterer (e.g. long-range Coulomb scattering, 
short-range defect scattering, etc.) under consideration 
with each scattering mechanism being independent. 

The finite-temperature conductivity is given by an ap- 
propriate thermal energy averaging within the Boltz- 
mann theory once r(e) has been calculated. The zero- 
temperature result is simply given by: 



-D(e F )T(e F ) 



(21) 



where the graphene Fermi velocity v F is assumed to be 
a constant (independent of momentum and density) and 
e F , the Fermi energy, is the chemical potential at T = 
0. The finite temperature chemical potential, fj,(n,T), 
is calculated self-consistcntly^ 5 - so that the net carrier 
density (induced by doping or an external gate) is n, and 
the gaplcssness of graphene automatically ensures that 
this procedure incorporates the thermally excited carriers 
(i.e., the only carriers present for intrinsic graphene as 
considered in Sec. II at the Dirac point with e F = 0) 
along with the induced carriers of density n. In this paper 
our main interest is the low carrier density regime where 
n is small so that the Dirac point behavior is accessed. 



A. Short-range disorder 

For short-range (or more appropriately, zero-range) 
delta scattcrcrs, we have 



\{V kk >)\ 2 =V?(l + cos0)/2, 



(22) 



where V& is the strength of the short-range disorder and 
the (1 + cos#)/2 factor arises from the matrix elements 
effect due to the pscudospin chirality of graphene (this 
chirality factor leads to the famous suppression of back 
scattering in graphene and also in surface states of topo- 
logical insulators). 

It is easy to show that short-range disorder leads to 
a carrier density independent conductivity a(n) oc V^~ 2 , 
and to an exponentially suppressed temperature depen- 
dence at low temperatures. In the high temperature 
limit, the resistivity due to short-range disorder increases 
by a factor of 2 compared with the T — valued 5 , 
and thus, short-range disorder by itself introduces weak 
metallic behavior in graphene with little temperature de- 
pendence at low (T -c Tp = ep/ks) temperatures and 
increasing resistivity at high temperatures (T — > oo). 
This is the same as what happens to just the Dirac point 
conductivity as discussed in section II. 



B. Long-range Coulomb disorder 

Unintentional charged impurity centers in the environ- 
ment are a major source of disorder for graphene on sub- 
strates. Although they are substantially removed in an- 
nealed SG samples (leading to the very high observed SG 
mobility), there are still some remnant random charged 
impurity centers on the SG surface which contribute to 
carrier scattering. For Coulomb disorder we have 



Vc(q) 



1 + cos ( 



(23) 



where V c (q) = 2ire 2 /nq 1 with k (= 1 for suspended 
graphene) as the background dielectric constant, is the 
2D Coulomb interaction and e(q) is the wave vector de- 
pendent static dielectric function of the free carriers in 
graphen o 10 ' 17 . 

The density dependence a(n) of conductivity due to 
Coulomb disorder is linear, a ~ n, and the preponder- 
ance of the observed linearity of c(n) on n is considered to 
be strong evidence for the importance of charged impu- 
rity scattering in determining graphene transport prop- 
erties. The temperature dependence due to Coulomb dis- 
order has been discussed elsewhere^ 5 -, and here we sum- 
marize the main findings for the discussion of our results 
presented in the rest of this paper. In the low temper- 
ature limit (T <C T F = ep/ks), one gets for Coulomb 
disorder 



a{T)=a [1- A 2 (T/T F ) 2 ] 



(24) 
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where A2 > 0. In the high temperature limit (T 3> Tp), 
which is the more appropriate regime for our considera- 
tion of transport near the Dirac point (n ~ 0), one gets 



o-(T) - B 2 (T/T F ) 2 , with B 2 > 0. 



(25) 



Thus, Coulomb disorder by itself predicts weak metal- 
lic behavior for T <C Tp and strong insulating behavior 
for T 3> Tp (which is the appropriate limit for the low 
density Dirac point regime). 



C. Acoustic phonon scattering 

In addition to the short-range and long-range disorder, 
which affect the SG conductivity at all temperatures (but 
with distinct density and temperature dependence in dif- 
ferent regimes), we also include resistive scattering by 
graphene acoustic phonons through the deformation po- 
tential coupling, which is operational primarily at higher 
temperatures (except at the Dirac point where it is op- 
erational at all temperatures). We note that the defor- 
mation potential coupling is rather weak in graphene, 
and therefore, the primary (essentially, the only) effect 
of phonon scattering is to introduce a weak metallic tem- 
perature dependence with the phonon-induced resistivity 
Pph ~ T at higher temperatures T > Tbg ~ 2fm p hkp, 
where v p h is the phonon velocity (i.e. speed of sound). 
Since kp oc y/n, phonon effects could affect the net SG re- 
sistivity at fairly low temperatures for low carrier density 
systems of our interest in this work. Since kp effectively 
vanishes at the Dirac point, acoustic phonons are opera- 
tional even at arbitrarily low temperatures near the Dirac 
point as Tbg tends toward zero. 

Since phonon scattering has already been considered 
in details theoretically elsewhere^ we show below the 
relevant "high-temperature" relaxation time for the de- 
formation potential coupling: 



1 



1 



D 2 



' J F Pn 



-{k B T), 



(26) 



"ph 



where D, p m are respectively the deformation potential 
coupling and graphene mass density. At very low tem- 
peratures, the phonon-induced relaxation time enters the 
Bloch-Gruneisen regime where p p h ~ T A and is negligibly 
smalU^. In our numerical results presented later in this 
paper, we use the full numerical solution of the Boltz- 
mann theory for calculating the phonon-induced resis- 
tivity which always becomes important above (a density 
dependent) characteristic temperature. 



D. Asymptotic behavior of SG conductivity 

We now combine the contribution to cr(T, n) = 
[p(T, n)] _1 from the three distinct scattering mechanisms 
described in subsections A, B, C above to discuss the 



asymptotic density and temperature dependence of the 
SG conductivity near the Dirac point. 

First we establish the counter-intuitive result that the 
conductivity around the Dirac point is always affected 
by phonon scattering even at arbitrarily low tempera- 
tures. Writing the effective carrier density n(T) around 
the Dirac point as 



n(T)« 7i + AT 2 , 



(27) 



where no oc V g is the gate induced extrinsic carrier den- 
sity and AT 2 (3> no) is the intrinsic Dirac point ther- 
mally excited carrier density (see Sec. II), we can define 
an effective Fermi wave vector: 



k F = v/7r(n + AT 2 



(28) 



where A = (ir/6)(kg /j) 2 are known T-indcpendent con- 
stants. Then, the Bloch-Grimeisen temperature Tbg 
above which phonon scattering effects are important is 
given by 



T BG = uJ P h(2k F )/k B = 2{v ph /v F )^/&T = dT, (29) 



where d = 2v p h/vpy / -K 2 (We note that Tbg is de- 
fined by the phonons with effective wave vector of 2kp 
since 2kp typically is the most resistive scattering process 
across the Fermi surface.) Now the condition for acous- 
tic phonons to contribute appreciably to the resistivity 
(e.g., p cx T) is that T > Tbg, implying d > 1. If d > 1, 
then the Dirac point resistivity remains unaffected by 
phonons to arbitrarily high temperatures, whereas d < 1 
implies that phonons contribute a linear resistivity down 
to low temperatures. It is easy to check that for actual 
graphene parameters, we find that d ~ 1CP 3 and thus 
d -C 1 is satisfied, implying that <td(T) is affected, in 
principle, by phonon scattering at all temperatures. We 
can estimate the crossover temperature scale T c for the 
low-density SG system to go from being 'insulating-like' 
dominated by Coulomb disorder to being 'metallic-like' 
dominated by phonon scattering to be T c ~ 2/(A c B p ) 1 / 2 
where A c is the coefficient for the T 2 dependence due to 
Coulomb disorder in Eq. (|T3|) and B p is the coefficient 
of the linear T-inverse term due to phonon scattering 



in Eq. II18]). It is easy to show that T c 



1/3 



and 



thus the crossover temperature increases with increasing 
Coulomb disorder in the system. For very pure SG sam- 
ples, T c could be very low, and it is in principle possible 
for the low-temperature Dirac point conductivity to show 
a transition from being insulating-like to being metallic- 
like as a function of decreasing disorder (i.e. n,), but 
this is by no means a localization transition - it is simply 
a crossover behavior driven by the competition between 
charged impurities and phonons. In general, the Dirac 
point conductivity would show complex non-monotonic 
temperature-dependent conductivity as is obvious from 
this analysis and from Fig. [T] 

For finite doping when no 3> AT 2 , the above argu- 
ment does not hold, and phonon effects on conductiv- 
ity are pushed to much higher temperature while at the 
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same time the temperature scale for the insulating behav- 
ior is substantially suppressed since Tp is now large [see 
Eqs. (|24|) and (f25j) above]. Then, the high-temperature 
behavior of <r(T) must always reflect the weak metallic 
(da/dT < 0) conductivity in the p p h oc T regime for 
T S> Tbg- Thus, at high carrier density both the in- 
sulating and the metallic temperature dependences are 
strongly suppressed leading to very weak temperature de- 
pendence of graphene conductivity as is well-established 
experimentally. How high in temperature must one go to 
manifest the weak metallic phonon-induced conductiv- 
ity obviously depends on the gate induced carrier den- 
sity. First, Tbg increases with increasing carrier density 
n since Tbg oc \/n. Thus, doped SG with high carrier 
density should reflect very weak temperature dependence 
except for the weak phonon-induced metallic behavior at 
high temperatures whereas the Dirac point conductivity 
(or more generally, low-temperature conductivity) should 
reflect strong (and in principle, nonmonotonic) temper- 
ature dependence in the conductivity. An observation 
of any strong temperature dependence in graphene (in- 
sulating, metallic, or nonmonotonic) therefore indicates 
the Dirac point behavior. 

For short-range disorder, as discussed in IIIA above, 
<t(T) has a weak T-dcpendence at low temperatures and 
metallic temperature dependence at high temperatures, 
and thus adding phonon scattering does not change the 
picture qualitatively. Thus, pure short-range disorder 
by itself can only introduce weak metallic temperature 
dependence in graphene transport properties with cr(T) 
decreasing with increasing T at high temperatures as 
phonons start playing a role. In addition, short-range dis- 
order does not manifest any density dependence in cr(n), 
and therefore, cr(T, n) would have little dependence on 
density and temperature (except at high temperatures) 
if the dominant resistive scattering mechanism is short- 
range defect scattering in conflict with all existing exper- 
iments. 

Long-range disorder, however, must dominate low- 
density transport since p cx 1/n for long-range disorder, 
and therefore, the Dirac point conductivity is necessar- 
ily limited by long-range disorder, which, as discussed in 
Sec. IIIB above, leads to a non-monotonic temperature 
dependence of weak metallic behavior for higher values of 
temperature. At low carrier densities (as well as in the 
presence of any remnant puddles) the low-temperature 
weak metallicity may be ignored, and the net tempera- 
ture dependence arising from long-range disorder (IIIB) 
and phonon scattering (IIIC) can be combined to give 

p(T) = (a + A.T 2 )- 1 + B p T, (30) 

where <jq subsumes the weak metallic temperature depen- 
dence at low temperatures and A c , B p (as defined above) 
depend respectively on the long-range Coulomb scatter- 
ing and acoustic phonon scattering as discussed above. 
For small ao, at the Dirac point, such a T-dependcncc 
leads immediately to a universal crossover temperature 



scale T c given by 

T c » {AB)~^ oc n\ ,3 *J*vy\y*D-W, (3 i) 

with T c being the characteristic temperature defining 
the crossover from the temperature power-law insulat- 
ing temperature dependence induced by Coulomb scat- 
tering to the higher-temperature phonon-induced weak 
metallic temperature dependence. We note that T c in- 
creases weakly with disorder (n,j) and r s , but decreases 
with increasing deformation potential coupling. Cleaner 
SG systems would thus manifest stronger phonon effects. 
But for high densities, when ao typically is large, and the 
phonon effects can only show up at very high tempera- 
tures T 3> Tbg(^ n o) w hh Tbg also being large, Eq. (|30|) 
immediately implies very little temperature dependence, 
except for p(T) — T for T > Tbg- Thus, away from the 
Dirac point, whence both Tp and Tbg & re large, SG con- 
ductivity should manifest weak temperature dependence 
whereas the observation of a strong temperature depen- 
dent conductivity is evidence for approaching the Dirac 
point in the system. 

In the next section, we provide our calculated numeri- 
cal results for SG transport properties using the full nu- 
merical solutions of the Boltzmann transport theory in- 
cluding long-range and short-range disorder and acoustic 
phonon scattering. 

IV. NUMERICAL TRANSPORT RESULTS FOR 
SUSPENDED GRAPHENE 

We consider three different scattering mechanisms in 
calculating the density and temperature dependent SG 
conductivity tr(n,T), or equivalcntly the resistivity p = 
1/cr: long range Coulomb disorder (ni), short range dis- 
order (nsVfi, where n$ is the short-range impurity den- 
sity and Vs is the strength of the short-range disorder) 
and acoustic phonon scattering (D). We assume D = 19 
eV throughout and assume that the long range and the 
short range disorder can both be taken to arise from ran- 
dom quenched point impurity centers located on the sus- 
pended graphene layer. 

In addition to conductivity (or resistivity) we also 
present results for the mobility p and the mean free path 
/ since these are quantities of considerable experimen- 
tal interest. In particular, the mean free path is often 
used by the experimentalists to operationally determine 
whether transport is ballistic or not - if I > L (where L 
is the system size) one nominally has ballistic transport 
(and our theory becomes inapplicable). Similarly, mobil- 
ity is an important physical quantity pertaining to the 
sample quality - typically SG samples should have high 
mobility because the amount of disorder is suppressed. 

In calculating the conductivity of extrinsic suspended 
graphene in the presence of finite doping (or gate induced 
carriers with a finite Fermi energy ep) we first generalize 
the theory of Sec. II to the finite doping case as discussed 
below (with no being the doping density). 
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FIG. 2. (color online) Temperature dependent electron den- 
sity n(T) [Eq. fl38j] and energy e{T) [Eq. (J36J)] as a function 
of temperature, T/Tf- 



The current density in the presence of an applied elec- 
tric field (E x ) is given by 



J x — E x 



where D(e) = ge/[2Tr(hv F ) 2 ] is the density of states of 
graphene with energy e = Hvpk and f(e) is the Fermi 
distribution function. Thus the conductivity becomes 



(33) 



To find a direct analogy of the conductivity with the 
parabolic dispersion, a = ne 2 (r) /m, we rewrite Eq. (|33[) 
as 



a(T) 



e 2 (r 



9<T) 
4irh 2 ' 



where 



JD(s)t(s) (-gg)& 

SD{e) (mf\ de 



(34) 



(35) 



which is exactly the same definition of the average scat- 
tering time for 2D parabolic band systems and e(T) is 
given by 

e(T) = J f(e)ds = /xo(T) + 1 In [l + e"^"] , (36) 

where Ho(T) is the chemical potential and at T = 
/Jo = £f- With the 2D parabolic energy dispersion 



e = (fikp) 2 /2m Eq. ([34|) becomes the 2D conductivity 
formula, er = ne 2 (r) jra. 

With the classical average velocity (v x ) the current 
density is given by 



J x = n(T)e(v x ), 



(37) 



where we use the total electron density at finite T instead 
of the zero temperature density tiq, and n(T) is given by 



(38) 



n(T) = J D(e)f(e)de. 
From Eqs. (JSHJ) and ([37]) we have 



en{T) 



E x 



(39) 



Then the mobility can be defined by (we note that we 
have used the standard notation /j to imply both mobil- 
ity and chemical potential which should not cause any 
confusion since they do not arise in the same equation in 
the text and it should be clear from the context whether 
mobility or chemical potential is being discussed) 



E x 



en(T) ' 



(40) 



Now we define the mean free path from the average 
scattering time as 



urr\ I \ h 2 hl ' F 
l(T) = v f (t) =<t—- 



e 2 g, 



(41) 



In mobility [Eq. (gDJ] and mean free path [Eq. (gTJ)] the 
total density and thermal energy at finite temperature 
are used instead of no and £p, i.e., 



and 



1{T) 



£(T)_ 

en 



h 2 Hvf 
e 2 g ef 



(42) 



(43) 



We note as mentioned already that there are two pos- 
sible alternative definition above for the mobility fi and 
mean free path I , depending on whether one uses the gate 
induced doping density no or the full carrier density n(T) 
including the thermally excited carriers. For T < Tp 
where Tp (= ep/ks) is always defined with respect to 
the T = carrier density induced by the gate, the two 
definitions are equivalent since n(T) ~ no = n. But, 
at very high temperatures (or very low doping density), 
n(T) 3> hq because T S> Tp. For ordinary graphene 
on substrates, where the puddle-induced density inho- 
mogeneity introduces a cut-off density of n c ~ 10 12 
cm -2 with a corresponding Tp ~ 1250A", by definition 
n(T) = riQ = n, and there is not much of a difference be- 
tween the two different ways of defining the mobility and 
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FIG. 3. (color online) Conductivity of suspended graphene 
corresponding to the experimental data of Bolotin et a/.—, (a) 
The calculated conductivity as a function of density for differ- 
ent temperatures T = 0, 100, 200, 300K (top to bottom) with 
m = 0.85 x 10 10 cm 2 and n s Vg = 1.5 (eVA) 2 . In (b) and (c) 
the mobility and mean free path are shown, respectively. The 
solid lines are calculated with temperature dependent n(T) 
and e(T), and the dashed lines are calculated with a zero tem- 
perature density no and the energy Ef- In (d), (e), and (f) 
we show <t, (i, and I at low densities (down to n = 10 7 cm -2 ), 
respectively. Note that as n — >• or T/Tf oo for a fixed 
temperature, u(T) oc <r(T)/T 2 and l(T) oc cr(T)/T. Thus, 
both fj,(n — > 0) and l(n — > 0) saturate at a finite temperature. 



the mean free path. But, for suspended graphene, where 
n c is small and thus very small values of rig are mean- 
ingful, the two distinct ways of defining fi and I make a 
big difference for small external doping density (no ~ 0), 
i.e., near Dirac point. 

To make the above point more explicit, we show in 
Fig. [5] our calculated n = n(T) and e(T) compared with 
no and ef% respectively. It is clear that for T/Tp > 1, 
there is a substantial difference between n(T) and no = 
n(T = 0). Since T F ~ 12K for n ~ 10 8 cm" 2 , close to 
the Dirac point, the two quantities (n and no) may differ 



a lot even at moderate temperatures. 
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FIG. 4. (color online) Conductivity corresponding to the ex- 
perimental data of Du et al£. (a) The calculated conductivity 
as a function of density for different temperatures T — 0, 100, 
200, 300K (top to bottom) with m = 5.0 x 10 10 cm 2 and 
nsVg = 14.7 (eVA) 2 . In (b) and (c) mobility and mean free 
path are shown as a function of density, respectively. The 
solid lines are calculated with temperature dependent n(T) 
and e(T), and the dashed lines are calculated with a zero 
temperature density no and an energy Ef- In (d), (e), and 
(f) we show a, fi, and I at low densities respectively. 

In Figs. I5HTU1 we show our calculated SG transport 
properties as functions of density and temperature ne- 
glecting all effects of density inhomogeneity or puddles 
- our theory should therefore be cut off at some very 
low doping density (< 10 9 cm -2 ) where puddles become 
relevant in high-quality SG. In Figs. [3]— [8j we show den- 
sity dependence for a few representative temperatures 
whereas in Figs. [5] and [TU] we show the calculated tem- 
perature dependence for a few fixed doping densities. We 
include phonon effects only in Fig. [9] and [10] since our 
main interest is low-temperature transport. In obtaining 
our numerical results, we focus on three published exper- 
imental SG works in the literature: Bolotin et al£-&, Du 
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FIG. 5. (color online) Conductivity corresponding to the 
experimental data of Mayorov et a/.—, (a) The calculated con- 
ductivity as a function of density for different temperatures 
T — 0, 100, 200, 300K (top to bottom) with m = 0.3 x 10 10 
cm 2 and nsVg = 1.5 (eVA) 2 . In (b) and (c) mobility and 
mean free path are shown as a function of density, respec- 
tively. The solid lines are calculated with temperature de- 
pendent n(T) and e(T), and the dashed lines are calculated 
with a zero temperature density no and an energy Ef- In (d), 
(e), and (f) we show a, a, and I at low densities, respectively. 



et al*^, and Mayorov et al^. Our goal is not data fit- 
ting or getting precise agreement with the experimental 
results since the details of disorder are unknown in the 
experimental systems, and we cannot rule out the pos- 
sibility of additional scattering mechanisms not included 
in our model being operational in the experimental sys- 
tems. What we are interested in is obtaining the broad 
qualitative features of SG transport data in our theory 
in order to critically assess the issues of ballistic versus 
diffusive SG transport and the access to the Dirac point 
at low carrier density. We mention that the authors in 
all three of these works^— interpret their data in terms of 
ballistic transport mainly by comparing their extracted 
mean free paths with the sample size. We will critically 
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FIG. 6. (color online) Calculated conductivity as a func- 
tion of density for different temperatures: (a) for long-range 
Coulomb potential with m = 10 10 cm~ 2 and (b) for neu- 
tral short range potential with nsVg = 5(eVA) 2 . Solid 
lines indicate a vs. n(T), and dashed lines indicate a vs. 
no = n(T — 0) or density induced by only gate voltage. 



examine the nature (ballistic or diffusive) of transport in 
these experiments. 

For each experiment, we choose a set of disorder pa- 
rameters as shown below based on the best overall semi- 
quantitative and qualitative agreement with the data, 
keeping these disorder parameters fixed for all the pre- 
sented results. The acoustic phonon scattering param- 
eters are standard and are taken to be fl = 19 eV; 
p m = 7.6 x 10~ 8 g/cm 2 ; v ph = 2x 10 6 cm/s. We use 
the following disorder parameters for each experiment. 



to. 



Du et al. rii = 6.5 x 10 



Bolotin et al. m = 1.2 x 10 



n s Vg = 15.7(eFA) 2 



10 cm- 2 , 



n s Vi = 1.5(eVAf 



Mayorov et al. n< = 0.3x 10 10 cm" 2 ; n 5 Vg = 1.5(e^A) 2 . 

We note that consistent with the experimental SG sam- 
ple quality, our disorder is the strongest (weakest) in 
Du (Mayorov) with Bolotin disorder being intermediate. 
This is consistent with the claimed high-density mobil- 
ity being - 5000, 000 cm 2 /Vs, - 200, 000 cm 2 /Vs, and 
~ 100, 000 cm 2 /Vs, respectively, in the three experiments 
(although the precise value of the sample mobility may 
not be a meaningful quality since the mobility depends 
on both density and temperature). 

In Fig. [3] — El we show our calculated a{n), fJ,(ri), and 
l{n) as a function of doping density n (alluded to hq 
above) for T = 0, 100, 200, 300K for the three experi- 
mental SG samples, respectively - we emphasize that for 
small values of no, where T/Tp > 1 condition may ap- 
ply, the alternative definitions for the chemical potential 
and the mean free path would lead to large quantitative 
differences since, as is obvious from Fig. [5J n(T) 3> n-Q in 
this regime. 

Our calculated cr(n) results for the three experimen- 
tal samples in Figs. [3H5] manifest similar qualitative be- 
havior with large quantitative differences because of the 
differences in the details of the underlying disorder. In 
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FIG. 7. (color online) (a) and (b) Calculated mobility as 
a function of density for different temperatures and for long- 
range Coulomb potential with m = 10 10 cm~ 2 . Here no in- 
dicates the density induced by the gate voltage and n(T) 
indicates the total density, i.e. the density from gate plus 
the density from thermal excitations. Solid lines represent 
Eq. ((40)) with n(T) and dashed lines represent Eq. ((42)) with 
no = n(T = 0) or density induced by only gate voltage. In 
(c) and (d) the mean free path are shown as a function of 
density. Solid (dashed) lines indicate Eq. ((41]) [Eq. (|43)) ]. 



particular, the following salient features of the results 
are consistent with the experimental findings in high- 
quality SG samples: (1) cr(n) manifests sublincar den- 
sity dependence, simulating a ~ y^n, over an extended 
density range, thus calling into question the experimen- 
tal interpretation of SG transport being ballistic^ based 
entirely on this sublinear density dependence; (2) at the 
lowest density, cr(n) is always limited by the long-range 
Coulomb scattering (with p oc n), but the competition 
between long-range and short-range disorder (which leads 
to the effective sublinear density dependence over an 
extended density range) and the existence of the low- 
density puddle-dependent cut-off (not included in the 
current theory) may mask this linear density dependence 
in high-quality SG samples where random charged im- 
purity disorder is presumably rather low; (3) at the low- 
est density, the system always would manifest insulating 
temperature dependence because of the dominance of the 
thermal excitation in the gapless system [this is obvious 
in the panel (d) of Figs.[3HS] near the Dirac point - there 
is a density-dependent crossover to the metallic behavior 
at higher carrier densities, emphasizing that the char- 
acteristic Dirac point insulating transport behavior is a 
high-temperature crossover behavior (which may not be 
apparent for T/Tp <C 1; (4) the calculated mobility ap- 



FIG. 8. (color online) (a) and (b) Calculated mobility as a 
function of density for different temperatures and for short- 
range neutral potential with ngVg = 5(eVA) 2 . Solid lines rep- 
resent Eq. (|40() with n(T) and dashed lines represent Eq. f|42[) 
with no = n(T = 0) or density induced by only gate voltage. 
In (c) and (d) the mean free path are shown as a function of 
density. Solid (dashed) lines indicate Eq. ((41]) [Eq. ((43]) ] . 



proaches ~ 5 x 10 5 , 10 5 , and 10 cm 2 /Vs respectively 
in the Bolotin, Du, and Mayorov samples close to the 
Dirac point, showing the unprecedentedly high qualities 
of these SG systems; (5) it is misleading to characterize 
the mobility (or the mean free path) using the gate in- 
duced density since this would produce erroneously large 
mobility and mean free path at low gate voltages, and 
in fact, would imply a divergent mobility (or mean free 
path) at the Dirac point - when the full density n(T) is 
used in defining the mobility (or mean free path) , the low- 
density mobility and mean free path saturate providing 
the correct characterization; (6) both definitions [using 
no or n(T)\ give identical mobility and mean free path 
values for high carrier densities (> 10 12 cm -2 ), as one 
expects since n(T) ~ tiq for high densities; (7) although 
broadly in qualitative agreement with the experimental 
data, there are important discrepancies between our the- 
ory and experiment in the details, most likely because 
of our neglect of other possible scattering mechanisms in 
the experimental systems; (8) the appropriate mean free 
path (at low densities near the Dirac point) varies be- 
tween ~ lOOnm (Du sample) and ~ lOOOnm (Mayorov 
sample), and therefore true ballistic transport measure- 
ments would require sample size < 0.1/im, and one must 
observe the sample length dependent conductivity to val- 
idate any ballistic transport behavior. 

Since long-range and short-range disorders affect 
cr(n, T) qualitatively differently (and the nature of dis- 
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FIG. 9. (color online) Temperature dependent conductiv- 
ity of suspended graphene corresponding to the experimental 
data of (a) Du et aL— and (b) Bolotin et aZ.— . The same pa- 
rameters used in Figs. [3] and [4] are used in this calculation. 
The solid (dashed) lines indicate the results with (without) 
phonon scattering. 



order in the experimental samples is not known based on 
any independent measurements), we depict in Figs. Eh 
IH1 the distinct theoretical dependence of conductivity on 
the long-range and short-range disorder separately (in 
contrast to Figs. [3H5] where both are included together 
in the theory) on the density for the SG system. To 
bring out the qualitatively different dependence of con- 
ductivity, mobility, and mean free path on no (the doping 
density) or n(T) at low density, we show in Figs. |6H1] de- 
pendence on both no and n(T) separately. In Fig. [51 
we show the calculated conductivity whereas in Figs. 
and [H] we show the mobility and the mean free path. In 
Fig. [6] we show u(n) for just long-range disorder or just 
short-range disorder with the two different density de- 
pendences [no or n(T)] showing quantitative differences 
only at low values of n (or equivalently, high values of 
T) with the two being identical (by definition) at T = 
since n(T = 0) = no- In Figs. [7J and HI we show the 
calculated mobility and mean free path for long-range 
(Fig. [7]) and short-range (Fig. [5J disorder with each case 
also providing the dependence on no and n(T). The im- 
portant qualitative conclusion from Figs. [BH5]is that one 
should always extract mobility and mean free path us- 
ing the correct total density n(T) rather than just the 
doping density no, particularly at low carrier densities 
because the extracted mobility and mean free path for 
the two definitions differ qualitatively as the Dirac point 
is approached with the distinction between the two def- 
initions being much larger for long-range disorder. Our 
work establishes that derived quantities such as mean free 
path and mobility, which involve an effective division of 
the experimentally measured conductivity by a density, 
are not meaningful for graphene (particularly at low den- 
sities, approaching the Dirac point) because n(T)/riQ di- 
verges at the Dirac point. This is not a serious problem 
for graphene on substrates because the puddle-induced 
cut-off density n c ensures that n(T) « no ~ n, but in 
high quality suspended graphene mobility and mean free 
path are meaningful only if they are extracted at high 
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FIG. 10. (color online) Temperature dependent conductiv- 
ity of suspended graphene corresponding to the experimental 
data of (a)(b) Bolotin et a/.— and (c)(d) Du et a/.—. The same 
parameters used in Figs. and 2] are used in this calculation. 
In (a)(c) the solid (dashed) lines indicate the results with- 
out (with) phonon scattering, and black (red) lines indicate 
the results for a density n = 10 10 cm -2 (n = 10 9 cm -2 ). In 
(b) (d) the individual contribution to the total conductivity 
(black line) is shown for a density n = 10 10 cm -2 . In (d) the 
non-monotonic behavior at high density does not appear due 
to the strong short-range potential scattering, but in high mo- 
bility samples (b) the non-monotonic behavior shows up due 
to the much weaker neutral impurity scatterings. 



density where n(T) « no. 

All the above results (Figs. ER3]) ignore phonon effects 
which are very weak in graphene and only affect high 
temperature transport. In Figs. 151 and [TU1 we show the 
explicit effects of acoustic phonon scattering in the the- 
ory by comparing results for a(T) including and exclud- 
ing phonons in the calculation at high (Fig. and low 
(Fig. [TO)) carrier density, and for the Bolotin et al££ and 
Du et alr^£ samples. In general, the phonon scattering 
effect is much stronger for the Bolotin et al. sample than 
the Du et al. sample because of the much higher quality 
(lower disorder) of the former. This finding is completely 
consistent with our theoretical analysis in Section II and 
III where we establish that the Dirac point conductivity 
would be affected by phonons even at rather low temper- 
atures for very clean samples with low values of n, . Our 
basic finding is that phonons introduce metallic temper- 
ature dependence at higher carrier density nullifying the 
intrinsic insulating temperature dependence arising from 
Coulomb disorder, but in general the insulating tempera- 
ture dependence remains quite strong upto the room tem- 
perature at low carrier density (Fig. 1 10[) in high-quality 
suspended graphene. We note that both Figs. [5] and [TU] 
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clearly show the very low-temperature (T/Tp <C 1) weak 
metallic T-dependence of cr(T) arising entirely from the 
Fermi surface effect which is more strongly manifested in 
the higher density (Fig. [5]) system. This again reinforces 
our claim that the insulating behavior of a(n,T), which 
is the hallmark of the Dirac point transport property, is 
much better studied as a high temperature phenomenon 
in low-density SG. This insulating behavior has clearly 
been observed by Bolotin et al, Du et al, and Mayorov 
et al, establishing that all three SG samples are reflect- 
ing intrinsic Dirac point transport behavior in their very 
high quality SG samples. Based on our results we con- 
tend that the observation of low density power-law in- 
sulating temperature dependence in graphene is a direct 
manifestation of the Dirac point behavior. 

V. CONCLUSION 

We have provided in this work a detailed theoretical 
study of the density and temperature dependent con- 
ductivity of low-disorder suspended graphene within the 
semiclassical Drudc-Boltzmann transport theory neglect- 
ing density inhomogeneity (i.e. puddle) effects. Our the- 
ory includes three independent scattering mechanisms: 
long-range Coulomb disorder, short-range (5-function dis- 
order, and acoustic phonon scattering. We establish, by 
comparing our detailed numerical results for the con- 
ductivity with three recent experimental studies^— , that 
the measured low-density conductivity in existing exper- 
iments on suspended graphene is approaching at least 
some aspects of the intrinsic Dirac point behavior. 

Some of our more important qualitative conclusions 
are: (1) the intrinsic Dirac point behavior is better man- 
ifested at higher (lower) temperatures (densities) staying 
above the puddle-induced characteristic density; (2) the 
observation of a power law insulating temperature de- 
pendence of conductivity is a direct manifestation of the 
Dirac point behavior; (3) at low doping densities, it is 
not meaningful to characterize the system using derived 



quantities (e.g. mobility or mean free path) because of 
the considerable ambiguity in which density (just the ex- 
trinsic doping density or the total density including ther- 
mal excitations) should be used in the definition of mo- 
bility (or mean free path); (4) the competition among 
long-range and short-range disorder plus phonon scatter- 
ing could lead to complex (and even non-monotonic) de- 
pendence of the conductivity on temperature and density, 
and it is not meaningful to conclude about the underly- 
ing nature of the transport behavior (ballistic or diffusive; 
localized or extended, etc.) based just on preconceived 
notions about the expected density and temperature de- 
pendence for various processes; (5) by improving sam- 
ple quality and reducing disorder, it should be possible 
to approach the Dirac point indefinitely through careful 
conductivity measurements in suspended graphene, pro- 
viding unique opportunities to study in the future many 
interesting effects not included in our theory (e.g, interac- 
tion, localization, ripple, flexural phonons); (6) phonons 
could affect the Dirac point conductivity in high-quality 
SG down to arbitrarily low temperatures since the Bloch- 
Griinciscn temperature becomes vanishingly small near 
the Dirac point - whether phonon effects will overcome 
the insulating temperature dependence due to Coulomb 
disorder depend on the details of the amount of disorder 
scattering effective in the system. 

We conclude by emphasizing that our results establish 
that how close in density one has approached the Dirac 
point can be estimated by seeing how high in tempera- 
ture the Coulomb disorder induced insulating tempera- 
ture dependence persists in a particular graphene sample 
(or paradoxically, how low in temperature the acoustic 
phonon effects persist if the graphene sample is devoid of 
Coulomb disorder causing the insulating behavior) . 
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